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I. ABSTRACT 

A model for a class of quasars and active jtalacttc nuclei la 
described In which a shock around a massive black hole randomizes 
the Infall kinetic energy of spherically accreting matter 
producing a non-theraal spectrum of high energy protons. These 
protons may be responsible for the secondary production (via t* 
decay) of the radio emitting high energy electrons and also of 
high energy y-rays (via s® decay and inverse Compton Interactions 
of the electrons). The correlation between radio and y-ray 
emission Implied by the model la In good agreement with obser- 
vations of 3C273. Observation of the flux of high energy neutrinos 
from quasars may provide a test for the model. 


II. INTRODUCTION 

Accretion onto compact oblects has long ^en considered as an 
energy source to power quasars (Lynden Bell 1*16*1) and Is thought 
to be the only mechanism capable of satisfying their enormous 
power requirements- Radio, X-ray and y-ray observations show that 
some quasars and active galactic nuclei emit roughly the same 
enersry per decade (l.e. per logarithmic energy Interval, d log-,g-) 
froa radio to X-ray energies and. In the case of 3C273. to y-ray 
energies. Relativistic particles required to account for the 
observed radio, X and y-ray emission may result from acceleration" 
in shocks via a first-order acceleration mechanism (Protheroe and 
Kazanas. 1*?S2). This siechanlsa, efficient (up to *18 percent; 
Axford IRRl) In transforming directed kinetic energy Into 




i . 

relativistic particle energy, has already been proposed to account 
for the observed cosmic ray spectncs (Bell 197fia, b: Axford, Leer 
and Skadron, 1977; Blandford and Ostriker, 1978; Cowslk and Lee, 

1981). For strong shocks, this process naturally produces a power 
law spectrum with index -2, l.e., equal energy per decade. The 
nuclear interaction of protons accelerated by this process results ^ 

in production of y-rays and secondary electrons responsible for 
synchrotron and inverse Compton emission. 


III. PARTICLE ACCELERATIOH AND RADIATION PROCESSES 

A shock around a massive black hole may serve both to 
randomize the Infall kinetic energy of the gas and also to produce 
a non-thenaal distribution of high energy (relativistic) protons. 

Because of the high photon energy densities 200 erg cm~^ for 

3C273 If the source region has dimensions less than 10^^ cm) and 

the magnetic fields present, electrons will suffer severe Cotspton 

and synchrotron losses, effectively prohibiting their direct 

acceleration. We shall, therefore, consider here secondary 

electrons resulting from nuclear collisions of the non-thermal 

proton distribution via plon production. Relativistic protons 

with kinetic energies much greater than the gravitational i 

potential of the black hole, do not readily fall into it. They 

can provide sufficient pressure to satisfy the lump conditions 

across the shock, the position of the shock being then fixed by 

the balance between Inlectlon and reaoval of relativistic protons *. 

behind It* Gamma-ray production via x°decay and via inverse | 

Compton scattering of the relativistic electrons arc considered. 

The spectrum of protons accelerated by,, the shock has a power 
law momentim dependence, l.e., Q (p) • K p ‘ (protons CeV~^ c s”^). i 

For a strong shock (compresslon^ratlo F • 2. The relati- 

vistic protons are assumed to be trapped within a radius comparable 
to the shock radius by magnetic fields which are tied to the 
Infalllng matter. The matter density in this region is assumed to 
be sufficiently high that the relativistic protons are depleted | 

mainly by nuclear Interactions with consequent plon production. • 

♦ 

The production spectnaa of y~rays, (E), resulting from j 

decay (Stacker 1971) has been calculated and Is plotted In ! 

Figure 1(a) where It Is normalized to the proton production 
spectrum at high energies, KpE~~- The production spectrum of 
secondary electrons Qg(E) (e'*’e~) from x-o-e decay has been 
calculated and Is also plotted In Figure 1(a). Above ' 1 CeV, 

Qg{S) ■ " 0-15 Qp(E) for T =• 2: the remaining energy goes 

Into neutrino production. Q (E) • 0.7 q (E). To obtain the absolute 
normalization of Qp(E) ^or 3C273, we sRall consider Its observed 
radio-optical spectrum. If we attribute the flat (F - v ) part 
of the soectrum from the radio to the vlslhle to synchrotron 
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radiation, then the total electron luainoslty, 2.3E Q ° (E), la 
approximately 2(1 + a) /a times the lynchrotron luminosity, where a Is 
the ratio of energy densities In the magnetic and radiation fields. 

Is the rate of production of electrons at high energies 
(electrons erg~^ s~^) and Is equal to Qg(E) plus the production 
spectrum of electrons and positrons resulting from photon-photon 
pair production interactions. 



Fig. Ir (a) Production 
spectra of electrons 
(Qg). y-rays from 
decay (<Vy^ y-rays 
from Inverse Compton 
Interactions for o “ 1 
(Qjc) relative to the 
production rate of 
protons. Separate curves 
are given for Q-j(. 
corresponding to 
$ ■ 0 and 1. (b) Total 

p^^uctlon rate of y-rays 
for • 1 relative to 
the rISlo synchrotron 
luminosity for different 
values of a. 


Assuming the electrons are trapped %dthln the shock radios, 

R, and lose energy by Inverse Compton and synchrotron radiation, 
the ambient spectrum of electrons, Ng(E), Is determined by 

m 111 

where, dE/dt ■ - 3.97 x 10~^ (I + a) ^rad 

the radiation energy density, taken to be 

inverse Compton Interactions, the y-ray production rate Is, 

Q--(E) K (E’) { f 2 n (E) o(E,e ,E')cde} dE' (2) 

IC 2 

where y « E'/mc^, n .(c) Is the number density of photons In the 
radiation field (phScons ca~^ erg”^), and o(E,e,E') Is the 
differential cross section for an inverse Compton interaction In 
which an electron of energy E* Interacts with a photon of energy 
e, boosting Its energy to E. 

We have solved equation (1) for Che two ertrese cases where 
the source is optically Chin ^ 0) thick (i = 1) to 

y-ray pair production Interactions. (^„^ Is the probability that a 
y-ray will annihilate with an X-ray phdton producing an electron- 





positron pslr). The resulting inverse Coopton T~ray production 
spectra are plotted in ^t^ure 1(a) for these two cases and a « 1. 
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tv. CAMMA-RAT SPECTRUM OP 3C273 


The observed y-rav spectrum will be the total Y-ray 
production spectrum (Q^ + modified by photon-photon pair 

production InteractlonsT To^calculate this effect, we oust Vnow 
the radius of the oblect, R. Ve shall obtain R from the frequency 
below which the synchrotron radiation becomes optically thick to 
synchrotron self absorption. Solvfn* equation (1), the number 
density of relativistic electrons of energy E within the source is 
n (E) •> n E , where n ■ (4.** x 10^^ L ^)/(aL^ ^R). 

t^;is spectrum, the frequency below whlch^XpSchrotron self 
absorption becomes Important, , is given by Pachoiczyk (1970), 

Vj - 2Cj(SCg)^^^ (3) 

where R is the component of magnetic field perpendicular to the 
line of^sight (B^ • ,R<B>), S ■ 1.3 R is the average distance 
through the source, and are constants. The synchrotron self 
absorption does not result in a sharp spectral break in 3C273. 

This nay be understood in terms of non-uniform radial distribu- 
tions of electron density and magnetic field (Condon and Dressel, 
1973). V, appears. to be in the ranee 3x10^^ - 1*1^? ,5* 3C273, 

hence we Mnd Ra • 5x10^' - 10^^ cm and <B>a • 14 - 75 


- 14 - 75 


gauss. 


Ve can now use these estimates of the radios to obtain the 

photon density in the X-ray region from the observed X-ray flux 

data and hence the effects of pair production interactions on the 

Y-ray spectrum. The optical depth to y-rays of energy E is 

t(E) - S / n j^(e) 0 (E.e) dc where n i^(e) is the number density of 

photons of eRergy €, and o(E,e) is the dlfferental cross section 

for pair production interactions of y-rays of energy E with 

photons of energy' , (Heitler, 1954), We have used X-ray data of 

Worrall et al. (1979; and Prlmlnl et al. (1979) sad the assumed 

radius to obtain n .(e). 

pb 

The energy spectrum of y-rays escaping from Che source is 
Chen down by a factor. (I - e ‘)/t, from that at production. For 
3C273 this factor is plotted as a function of energy in Figure 2(a) 
for various assumed radii- We sec that the source is optically 
thick to y-rays above ' 100 MeV for radii in the range obtained 
earlier. We have therefore plotted in figure 1(b) the y-ray 
spectrum at production (Q Q,-) ^or ^ » i and for values of 

a ranging from I to 1/16. ^The observed y-?av lusinosity (Herasen 
et al.. 1981) is compared In Figure 2(b) with our prediction 
taking into account photon-photon absorption and is found 
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Pig. 2: (a) Redaction factor 

applied to predicted optically 
thin Y-ray spectra for 3C273 
for various assumed radii in 
order to obtain the expected 
Y»ray emission, (b) Camoa-ray 
Tuainosity of 3C273 based on 
the observed spectnsa (Rermaen 
et al., 1^81) compared with 
the luminosity predicted for 
a ■ 1/10 and for the range of 
radii obtained as described in 
the text. 


to be in good agreement vlth the prediction for a •• 0.1 and R " 
10^^ ca. This radius is consistent with that »e obtained earlier 
by considering the radio spectrum sosse 12 decades Tower In energy. 

V. COKCLUSIOHS 

^ have discussed a model for a class of quasars and active 
galactic nuclei, gauged to 3C273. It focuses on the high energy 
particles in these oblects which are responsible for the radio and 
Y-ray emission and predicts a correlation between these two 
extr^e parts of the speetrus which appears to be In good 
agrees^nt with observations of 3C273. The model provides 
estimates of the central s<nirce radius and Its magnetic field. If 
the total observed luminosity (L . “ 2x10^ erg s ) is due to 
gravitational energy release at che^surface of the shock, and the 
relativistic particle pressure balances the accretion ram 
pressure pv^, then the radius Inferred above Implies a mass 
of *■ 10^® Mq, and an accretion rate of '■lO^ yr“^ . 

The estimated radios is compatible with variability time 
scales of the order of several months which have been observed 
(Worrall et al. 1979). Variability reported on tsiuch shorter time 
scales {'■ I day) (Marshall. Warwick and Pounds, 1981) would 
clearly be Incompatible with our model. However, the statistical 
significance of these data Is low and they may also be compatible 
with a much longer time scale variability. It is important to 
point oat here that a radius of ^ 1 light day would render 3C273 
completely ooaque to Y“rays of energies greater than a few MeV, 
contrary to observation. 



Pair production by high energy y-taya and production of 
electrons by plon decay throughout the source voluae avoids the 
Coapton Catastrophe (Hoyle, Burbidge and Sargent, 1466) and allows 
a higher energy density in the radiation field than in the 
magnetic field (i.e. a < !)• The model also provides a situation 
sii^lar to that considered in the literature for the formation of 
lets (Smith et al. I960), i.e., a hot gas (the relativistic 
particles) underneath a cold gas (the infalling cold gas). 


Based on this e»del, one would expect emission of high energy 
Y-rays from other quasars and active galactic nuclei. So far, Cen 
A is the only other source detected in high energy Y~rays 
(Grlndlay et al., 1975). The observed flux is compatible with 
that predicted by a sitaple model for Cen A (Razanas and Protheroe 
1982) simitar to that discussed here for 3C273. 


Finally, the model predicts a higher flux of high energy 
neutrinos chan in models where electrons are directly accele- 
rated. ife predict a flux of 1.4 x 10 x (l+2a) neutricos (all 
types) ca~^ s“^ above I TeV for 3C273. This may be iri.thln the 
detection capability of the proposed DOMAHD neutrino telescope 
(Stenger, 1981) and thus asy provide a test of Che teodel. 
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